Possible explanation for the absence of bilayer splitting in YBCO 
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It has been claimed that the absence of the bilayer splitting in the high-temperature superconductor 
YBCO is a strong experimental indication that there are no coherent quasiparticles present in the 
CuO planesu. We study a pair of strongly correlated planes which are connected by a hopping 
transfer integral t± in the limit of large in-plane coordination number. The effect of the correlation 
is incorporated in a dynamical mean-field theory, where the Weiss field is determined by a two-site 
Hubbard Hamiltonian. We have solved this problem by numerical techniques and present results 
for the spectral function p for relevant parameters of the model. For small t± we find the coherent 
hoping between planes to be proportional to the in-plane quasiparticle renormalisation. 

PACS: 71.10.-w,71.10.Ay,71.10.Fd,71.28.+d,71.30.+h,74.25.Jb,74.72.-h 
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The transport properties of the high temperature su- 
perconductors show an unusual behavior in many re- 
spects. Andersontl has emphasized the difference between 
the transport properties in a CuO plane and perpendic- 
ular to the CuO plane stressing in particular the optical 
conductivity, in the c-direction perpendiculajj-to the CuO 
planesB. Also photoemission experimentsau in YBCO 
show only a small if any splitting of the bonding and 
antibonding CuO plane bands due to coherent hopping 
between the CuO planes. This is even more surprising 
in view of the relatively large value of the bare hopping 



matrix element, t±_ rs 0.3 eV compared to t 



Ml 



leV, 



estimated from bandstructure calculationso. The exper- 
imental data have been interpreted as a signal—for the 
absence of coherent quasiholes in the CuO planed. 

The question of how interactions affect the coherent 
transport between the planes has been studied mostly on 
one dimensional strongly correlated systems because in 
that case it is clear that the excitations above the ground 
state are spinons and holons as in the original proposal 
by Anderson. Two coupled chains (ladders) have been re- 
cently studied in detail Jpsfjaany-groups with numerical 
and analytic techniques&BtfflliJO. The result was that 
the hopping matrix element t± is a relevant perturba- 
tion for physically relevant values of the Luttinger liquid 
exponents and that the transport between the chains is 
coherent. Recently it has been claimed from an exact di- 
agonalization study that for longer range interactions in 
the chains the coherent hopping between the chains can 
be reduced to zeroJjor a finite value of the bare hopping 
matrix element t±tn. 

The problem of two coupled CuO planes is clearly 
more difficult since at this point the ground state and 
the excitation spectrum even for a single plane of the 



two dimensional Hubbard model away from half filling 
are not known. The problem has been studied with niA=. 
merical techniques like Quantum Monte Carlo (QMC)E^I 
and the fluctuation-, exchange approximation (FLEX) 
approximatiorotJ'tj. The subtle features of the dynam- 
ics of the quasiparticles in the coupled planes can cur- 
rently not be resolved in the QMC due to finite size-effects 
although the resultsJac pae plane look encouragingu. The 
FLEX calculationfOO'Ej are able to resolve the dynam- 
ics and yield surprisingly good results when compared to 
the QMC calculation but this approximation is uncon- 
trolled in the regime where the interaction is large. 

In this paper we present a calculation of the spectral 
function and the reduction of the coherent hopping in an 
approximation where the local correlations are taken into 
account exactly. We are in particular interested in the 
dynamics of the quasiparticles. We consider a model of 
two strongly correlated planes. We neglect the Coulomb 
interaction between the two planes. The motion of the 
holes in each plane is described by a Hubbard Hamilto- 
nian. 
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where t is the hopping matrix element between two near- 
est neighbors i and j, fi is the chemical potential, c\ i(T 
is the electron creation operator on site i with spin a in 
plane a (a = 1, 2) and n a ia is the electron number oper- 
ator for spin a at site i in plane a, n a i„ — <^ aia c aia . We 
denote the electron creation and annihilation operators 
in plane 1 by and c xi respectively and in plane 2 by 
and Cos- The two planes are coupled by the hopping 



1 



term 



H 



(2) 



with the hopping matrix element t± connecting site i in 
plane 1 with site i in plane 2. The complete Hamiltonian 
for the bilayer is then H = Hx+ H 2 + Hi 2 . An exact solu- 
tion of this problem at this point in time is clearly impos- 
sible since already Hi and H 2 cannot be treated exactly 
in two spatial dimensions. However some progress can be 
made by introducing the so called dynamical mean field 
approximation in which the local correlations a*e4ireated 
exactly. For a review of this method see e.gJlZHla. The 
important point is that the dynamics of the quasiparti- 
cles hopping onto a given site and leaving a given site 
is kept. To apply this approximation to the problem at 
hand the number of neighbors z of each site in a given 
plane (1, 2) to its neighbors in the same plane has to go to 
infinity, see Fig. |l|. In this approach the phase informa- 
tion of the electrons hopping from a site j to a given site 

1 in the same plane is lost because the phase information 
is averaged of over the large number of sites j connected 
to i. The hopping term between the planes has to be 
treated separately because there is only one site in plane 

2 connected to a given site in the plane 1. 




FIG. 1. Site of the two planes connected by nonzero hop- 
ping matrix elements 

The non interacting (U = 0) Green's function G^ 1 is 
defined as [G (fc,r)] aa , = - < T T c QfeQ (r) C t, fc/3 (0) > with 
a, a' G {1, 2}. It is given by 



G^(k,uj) = u-H(U = 0) = 



OJ - fife t± 

t± uj- e k 



(3) 



In the paramagnetic phase the self energy E(fc,u;) has 
only two independent components. We denote them by 
En and Ej_ in the present basis 



E|| E_l 

S_l En 



(4) 



The Hamiltonian is diagonal in the bonding - antibonding 
representation. The full Green's function can be obtained 
from the Dyson equation 



G^ 1 {k,uj) = G^fcu) - E(fc,w) 



(5) 



We are now in the position to derive the dynamical mean 
field equations for the coupled planes. In the limit of large 
coordination number the. self energy becomes local and 
therefore fc-independento. The only fc-dependence of the 
Green's function is left in the free electron dispersion 
efe. This fc-dependence can be absorbed in the density 
of states D(e) = J2k — £fe )> which of course depends 
on the type of lattice in which the electrons are moving. 
For the hypercubic lattice D(e) — exp(— (e / V2t) 2 ) / t^/2n . 
The full local Green's function is therefore given by: 



G(w)= / de k D(e k ) (Gq w) - 



(6) 



In the dynamical mean field approach C/q 1 (r) plays the 
role of a Weiss field. The main point of this approach is 
that the local self energy of the full problem is determined 
by the local impurity modeled: in our case a two impurity 
model. It has the action 

S = f dr f dr'Y^^{T)g^{r - r')^(r') 
Jo Jq a 

/ dTi, a1 (T)^ a1 (r)^ aL (T)^ al (T) (7) 



where ip a — (tpia-, V , 2<t) is a spinor of Grassmann vari- 
ables. We only consider the paramagnetic state. The 
Weiss field Qq 1 ^) is therefore diagonal in spin space. 
Using Eq. ^ we obtain the full local Green's function in 
the bonding, anti-bonding basis: 



G(w) = 



D + (u) 
£>_(w) 



(8) 



where D±(w) = D(uj - E||(w) ± (tj_ - E_l(u;))) is the 
Green's function of the anti-bonding and bonding band 
respectively and 



D(uj) = / de 



D(e) 



(9) 



is the Hilbert transform of the density of states D(e). 

We solve the two impurity model, pEq. @i using the 
"iterated perturbation theory" (IPT)u which has been 
shown to yield qualitative agreement with QMC and ex- 
act diagonalization approaches for interaction strength 
up to U/t ~ 4. The self energy is approximated by the 
first two terms of a perturbation expansion in the inter- 
action U. However these two terms are calculated with 
the local Green's function Qq which is then determined 
sclfconsistcntly. For doping levels away from half filling 
we used-,an extension of the IPT due to Kajueter and 
KotliaiM 
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In Fig. g the spectral functions of the bonding and an- 
tibonding band are shown as a function of frequency for 
a U = 4, a bilayer hopping of t± = 0.3 and a filling of 
n = 0.4375. All energies are measured in units of the 
in-plane hopping matrix element t. For each band the 
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FIG. 2. Spectral functions for the bonding (solid line) and 
anti-bonding (dashed line) band. The parameters are U = 4, 
n = 0.4375 and t± = 0.3. The dotted vertical line marks the 
Fermi energy. 

lower and upper Hubbard band can be clearly identified. 
The quasiparticle peak close to to = for each band is 
slightly shifted away from zero as expected from the non- 
interacting limit. The first interesting observation is that 
the spectral weight in the lower and the upper Hubbard 
band is not symmetric anymore even at half filling. The 
spectral density of the upper Hubbard band of the bond- 
ing band and the lower Hubbard band of the anti-bonding 
band are suppressed. With increasing hopping between 
the planes the quasiparticle peaks move further apart and 
the features associated with the Hubbard bands become 
less and less pronounced. For very large values of ij_ the 
nonintcracting density of states of the two bands is recov- 
ered. The small features in the spectral function which 
are associated with the formation of so called shadow 
bands as obtained in the FLEX calculation^ can not 
be found in our calculation, because the finite momen- 
tum antiferromagnetic fluctuations are missing from our 
calculation. However the redistribution of weight from 
the Hubbard bands to the quasiparticle peak and vice 
versa has not been resolved in the FLEX perhaps be- 
cause it is inherently a weak coupling approach. At half 
filling both bands are filled up to co = 0. The self energy 
of the bonding and antibonding propagator is given by 
E||(w) ± £j_(w). For small values of t± the quasiparticle 
weight is strongly reduced. 

Next we discuss the reduction of the coherent hopping 
between the planes. There are two measures of coher- 
ent hopping: one is AE, the splitting in frequency be- 
tween the quasiparticle peaks in the bonding and anti- 
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FIG. 3. The splitting AE/2t±_ (AK/2t±) of the coherent 
hopping between planes as a function of tj_. The parameters 
are (7 = 4 and n = 0.5. The dashed line refers to the splitting 
in energy, the solid line to the splitting in efe. The spectral 
weight of the quasiparticle peaks is shown as dashed-dotted 
line. 

bonding bands; the other is AK = 2(t± — T,±(oj = 0)), 
the splitting in momentum space between the bonding 
and antibonding bands at the Fermi surface. We con- 
sider AE/2t± first. Its t± dependence is shown in Fig. 0. 
In the noninteracting case AE assumes the value 2t±. 
For small t± the ratio AE/2tj_ is clearly reduced from 
the noninteracting value of 1. However unless one drives 
the system through the Mott transition at very large val- 
ues of U where the planes themselves are insulating the 
reduction factor is nonzero. In the t± — > limit the re- 
duction factor can be analytically related to the reduction 
of the bandwith in the plane. The off-diagonal part of 
the self-energy Ej_ is proportional to for small t±. The 
position of the quasiparticle-peak is given by the poles of 
D±(z) in the complex plane. For small t± and doping 
not too far from half filling the position of the pole is 
given by the zero of lu — Re £j|(w) ± t± + 0(t\) = 0. 
The selfenergy does not depend on the momentum. The 
reduction of the bandwidth is therefore solely given by 
(1 — <9E||(a;) / 'duo)^ 1 which is the same as the reduction 
of the bandwith due to interaction in a single plane. In- 
terestingly the coherent hopping is seen to be enhanced 
for larger values of t^. This effect can be understood in 
terms of a level repulsion of the two site Hubbard model 
in the regime where t± U t. For very large values 
of tj_ the splitting asymptotically approaches the value 
2t± as one would expect from the nonintcracting system. 

In Fig. [| we plotted the spectral density A(ek,uj) at 
zero temperature as function of frequency for parameter 
values ij_ = 0.3 and U — 4. The wavevector dependence 
is only implicitly contained via e^, because of the local 
nature of dynamical mean field theory. For close to 
zero we observe one sharp peak associated with the bond- 
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FIG. 4. Single particle spectral function A(ek,io) vs. u> for 
different values of ek- The parameters are t± = 0.3, U = 4 
and n = 0.4375. 



ing band just below the Fermi level as one would expect 
(The bonding band actually crosses the Fermi level at 
£k ~ 0.15). When going to lower energies this maximum 
moves downwards in the spectrum and very rapidly be- 
comes broader. The broadening is due to strong correla- 
tions which manifest themself in a rapid increase of the 
self energy away from u> = 0. For ek ~ —0.55 a new sharp 
peak due to the antibonding band appears. The peak 
of the bonding band, however, is already rather broad 
at this energy. The splitting AK between the ek val- 
ues where the bonding and antibonding bands cross the 
Fermi surface is thus AK w 0.7 so for these parameters 
AK/2t± w 1.1. 

In conclusion we have studied the problem of two cou- 
pled planes in the dynamical mean field theory in the 
paramagnetic state as a function of interaction, doping 
and interplane hopping matrix elements. Within our ap- 
proximation the coherent hopping between the planes is 
strongly reduced due to the interaction however it stays 
nonzero as long as the planes stay in the paramagnetic 
phase. It is still an open question if short range antifcr- 
romagnetic correlations which are completely neglected 
in our approach can suppress the coherence between cou- 
pled planes even further. For intermediate values of the 
interplane coupling the coherent hopping motion of the 
quasiparticles between planes is increased by the inter- 
action. The spectral weight in the Hubbard bands is 
redistributed in an unsymmetric way even at half filling. 
Our approach can be easily extended to include more 
realistic density of states. It is interesting to examine 
the stability of other groundstates like for example the 
antifcrromagnetic groundstate. 
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